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Momentum dependent local-ansatz (MLA) wavefunction describes accurately electron 
correlations from the weak to intermediate Coulomb interaction regimes. We point out 
that the MLA can describe the correlations from the weak to strong Coulomb interac- 
tion regimes by modifying the starting wavefunction from the Hartree-Fock (HF) type 
to an alloy-analogy (AA) type wavefunction. Numerical results based on the half-filled 
band Hubbard model on the hypercubic lattice in infinite dimensions show up that 
the new wavefunction yields the ground-state energy lower than the Gutzwiller wave- 
function (GW) in the whole Coulomb interaction regime. Calculated double occupation 
number is smaller than the result of the GW in the metallic regime, and is finite in the 
insulator regime. Furthermore, the momentum distribution shows a distinct momentum- 
dependence in both the metallic and insulator regions, which are qualitatively different 
CN ■ from those of the GW. 
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The variational theory is a useful tool for studying the ground-state properties of 
correlated electron system such as the magnetism, the heavy-electron behavior, the 
metal- insulator transition, and the high-temperature superconductivity. 1 " 16 ^ There the 
form of the wavefunction is crucial to describe electron correlations. The Gutzwiller 
wavefunction (GW) is one of the popular wavefunctions, because of its simplicity and 
applicability to realistic systems. 1-3 ) It describes electron correlations by making use of 
a projection operator IIj(l — gh^n^) onto the Hartree-Fock (HF) wavefunction \<f>o); 



* E-mail address: kl08609@eve.u-ryukyu.ac.jp 
^yok@sci.u- ryukyu.ac.jp 



J. Phys. Soc. Jpn. 



LETTERS 



I^gw) = [IL;(1 — g^i^i])] I0o)- Here hi a is the number operator for an electron on site 
i with spin a, variational parameter g reduces the amplitudes of doubly occupied states 
on local orbitals. The local-ansatz approach (LA) proposed by Stollhoff and Fulde is an 
alternative method which is simpler than the GW 4 " 7 ' The LA wavefunction takes into 
account the states expanded by the residual Coulomb interactions {Oi} = {^n^^n^}; 
I^la) = [ riiC 1 -VLA°i)] \<fio)- Here 5n ia = n ia - (n ia )o, (n ia ) being the HF average of 
electron number on site % with spin a. The amplitude 77 LA is determined variationally. 

Although the GW and the LA have been applied to a number of correlated electron 
systems, they are not sufficient for the description of correlations from the weak to 
strong interaction regimes. In fact, both the GW and the LA do not describe exactly 
the weakly correlated regime due to the use of limited Hilbert-space. Furthermore in 
the strong Coulomb interaction regime, the GW yields the Brinkman-Rice atom (i.e., 
no charge fluctuation on an atom) instead of the insulator solid in infinite dimensions 
(d = 00). 17 ) In order to overcome the difficulty in the weak Coulomb interaction regime 
and to improve the behaviors in the intermediate Coulomb interaction regime, we have 
recently proposed the momentum-dependent local ansatz wavefunction (ML A). 18 ' 19 ) 
The MLA wavefunction controls the amplitudes of momentum-dependent two particle 
states to be best, and much improves the Fermi liquid properties of the GW and the 
LA in those regimes. 

The wavefunction which accurately describes the strong Coulomb interaction regime 
has been proposed by Baeriswyl. 9-12 ) The Baeriswyl wavefunction (BW) is constructed 
by applying a hopping operator T onto the atomic wavefunction l^oo); I^bw) — 
e~ T)T \^ 00 ). Here T = — X)i j a^ij a la a j^ * s the kinetic energy operator, denotes the 
transfer integral between sites i and j, a ia (aj CT ) being the creation (annihilation) opera- 
tor for an electron on site % with spin a. The operator e~ vT with a variational parameter 
rj describes electron hopping from the atomic state and suppresses the configurations 
with high kinetic energy. Although the BW describes well the insulator state in the 
strong correlation regime, it is not easy to describe the metallic state from this view- 
point. There is no other variational wavefunction which is suitable for both the weak 
and strong Coulomb interaction regimes and can be analytically treated, as far as we 
know. 

In this letter, we point out that a modified MLA wavefunction which starts from the 
alloy-analogy (AA) wavefunction instead of the HF one describes the strongly correlated 
regime reasonably so that one can go beyond the GW in both the weak and strong 
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Coulomb interaction regimes. The wavefunction provides us with a new tool to describe 
various systems with intermediate Coulomb interaction strength such as the cuprates 
and the iron pnictides superconductors. 

We adopt in this letter the single-band Hubbard model as follows: 1 ' 20,21 ) 

H = ^(e - fijhia + ^ *ij" a !r a Jff + U ^2 ^T™4- (!) 
ia ijtr i 

Here eo (/x) is the atomic level (chemical potential), tij is the transfer integral between 
sites i and j. U is the intra-atomic Coulomb energy parameter. 

The construction of the MLA wavefunction is rather simple. We expand first the 
LA wavefunction up to the first order with respect to the Coulomb interaction strength 
U and observe that each coefficient of the two-particle excited states in the first-order 
wavefunction is momentum independent in contradiction to the exact result. Therefore, 
we introduce a new set of local operators {Oi} with the momentum-dependent varia- 
tional parameters rjkfak^ki which yields the correct weak Coulomb interaction limit: 

Oi= (fciKX^iX^KMfo) ^fc 2 fcife 1 ^(4^ a fc24)^(4it afe it) • ( 2 ) 

kik' 1 k2k' 2 

Here (i\k) = exp(—ik ■ R^/y/N is an overlap integral between the localized orbital 
and the Bloch state with momentum k, Ri denotes atomic position, and N is the 
number of sites. a\ a {a ka ) denotes a creation (annihilation) operator for an electron 
with momentum k and spin a, and 8{a\, a aka) = a k'a a k<r ~ ( a !'cr a fc°-)o- We then construct 
the MLA wavefunction as follows. 18,19 ^ 

|^MLA) = n( 1 -^)l^)- ( 3 ) 

i 

The best wavefunction is chosen by controlling the variational parameters in the mo- 
mentum space. Hereafter we refer the wavefunction (3) to the MLA-HF wavefunction. 

The MLA-HF wavefunction does suitably describe the electron correlations in the 
weak and intermediate Coulomb interaction, but it cannot suppress loss of Coulomb 
repulsion in the strong interaction regime. To improve the difficulty we propose here to 
change the starting wavefunction from the HF wavefunction to the alloy-analogy (AA) 
wavefunction which is suitable in the strong Coulomb interaction regime. 

The concept of the AA can be traced back to Hubbard's original work on electron 
correlations. 22 -* In the strong Coulomb interaction regime, electrons with spin a move 
slowly from site to site due to electron correlations, and therefore should feel a potential 
U instead of the HF average potential U (n i _ cr )o, when the opposite spin electron is on 
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the same site. Hubbard regarded this system as an alloy with different random potentials 
e + U and e . The AA Hamiltonian is then defined by 

H A A = ^(e - H + Uni_ a )h ia + ^ t ij a ia a j<r ~ U ^2 n it n H- ( 4 ) 
icr ija i 

Since the motion of electrons with opposite spin are treated to be static in the AA 
approximation, related operators {n ia } are regarded as a random static C number n ia 
(0 or 1). Each configuration {n ia } is considered as a snapshot in time development. 

We adopt the AA ground-state wavefunction 0aa for the Hamiltonian Haa, and 
propose a new ansatz, which we call the MLA-AA wavefunction, as follows. 

|^mla-aa) = Y[(l - OOI^aa). (5) 

i 

The local operators {Oi} with variational parameters r) K ' 2K2K > iKl have been modified as 
Oi= Yl ( /c il i )<*l«i><'4l i ><*l«2> ^/caiciKi^Ca^a^^^a^a^) . (6) 

Here a\ a and a Kcr are the creation and annihilation operators which diagonalize the 
Hamiltonian Haa, and ^{a^^a^) = a\, a a KCT — {a\, a a Ka )Q. It should be noted that the 
MLA-AA wavefunction reduces to the MLA-HF by replacing the random potential 
Urii- a with the HF one, i.e., U (0o|^i-o-|0o), so that \I/mla-aa and ^mla-hf are mutually 
connected to each other via a suitable parameter which interpolates between the two 
wavef unctions. 

We can obtain the ground-state energy for the MLA-AA wavefunction within the 
single-site approximation (SSA) taking the same steps as in the MLA-HF. 18, 19 ) The 
correlation energy per atom is then given by 

-(d t ] H) Q -{Hd l ) + (d l ] Hd l ) 

e c — _ i _ • {<) 

1 + (CVOi)o 

Here H = H — (H) . Note that the average (~) is now taken with respect to the AA 
wavefunction. The total energy per atom should be obtained by taking the configu- 



rational average; (H) = (H) + Ne^. Here the upper bar denotes the configurational 
average. Each term in the correlation energy can be calculated by making use of Wick's 
theorem. 

We obtain a self-consistent equation with variational parameters using the minimum 
energy condition. 
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- $^(! ~ /( g 4t))/(w)^«44« + ~ /( g 4t))( 1 ~ /(e«u))^ 4 «24fei = ^ • ( 8 ) 

Here A£^ K2K ' iKl is the two-particle excitation energy given by AE K > 2K2lt ' iK1 = e K ^ - 
+ e Kit — e Kit- e ro denotes one-electron eigen value energy for Haa, and e Kcr = 
e«;cr— A*- /(e) is the Fermi distribution function at zero temperature. To solve the equation 
approximately, we make use of an interpolate solution which is valid in both the weak 
Coulomb interaction limit and the atomic limit. 18 ) 

. Ufj . . 

Vh' 2 K2K , 1 k 1 \Vi e c) — ; ^ _ • 

Here 77 = [1 — 77(1 — 2 (71^)0) (1 — 2(714)0)]. The best values of fj and e c are determined 
again variationally. 

In the SSA, the correlation energy e c can be expressed by the local density of states 
Pi a (e) for the A A Hamiltonian (4). Because there is no translational symmetry due to 
random potential, we calculated p«j(e) by means of the Coherent Potential Approxima- 
tion (CPA). 23 ' 24 ) 

To examine the validity of the MLA-AA as well as MLA-HF, we have performed the 
numerical calculations for the half-filled band Hubbard model with nearest neighbour 
transfer integral on the hypercubic lattice in infinite dimensions, where the SSA works 
best. 25 ' 26 - 1 We assumed here the non-magnetic case. In this case, the density of states 
for non-interacting system is given by p(e) = (1/y/n) exp(— e 2 ). 25 ) The energy unit is 
chosen to be J dep(e)e 2 = 1/2. The characteristic band width W is given by W = 2 in 
this unit. 

Figure 1 shows the results of the ground-state energy vs Coulomb interaction energy 
curves. In the weak Coulomb interaction regime (U/W < 1), the total energy of the 
MLA-HF is lower than the GW. On the other hand, we observe that the MLA-AA gives 
lower energy in comparison with the GW in the strong Coulomb interaction regime 
(U/W > 1). We obtain the critical Coulomb interaction U c2 = 3.40 at which the 
effective mass diverges. But before U approaches U C 2 we find that the AA state showing 
the insulating state is stabilized, and the metal-insulator transition occurs at the critical 
Coulomb interaction U c = 3.26. The transition is the first order in the present approach, 
and is consistent with the result of numerical renormalization group method (NRG), 
although the calculated U c2 is somewhat smaller that obtained by NRG ( i.e., U c2 = 4.1 
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Fig. 1. The energy vs Coulomb interaction energy U curves in the MLA-AA (solid curve), MLA-HF 
(dashed curve), and GW (thin solid curve) for the electron number n = 1.0. 
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Fig. 2. The double occupation number (n^n{) vs Coulomb interaction energy U curve for the elec- 
tron number n = 1.0 in the MLA-HF (solid-dashed curve), MLA-AA (dashed-solid curve), and GW 
(thin solid curve). The arrow shows a jump from the metallic state to the insulator at U c = 3.26. 



). 27 ) The MLA-HF leads to the total energy lower than that of the GW up to the 
critical Coulomb interaction U c and the MLA-AA gives the same behavior in the range 
U > U c . More important is that the MLA scheme gives lower energy for overall Coulomb 
interaction regime and therefore can overcomes the limitation of the GW. 
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We present in Fig. 2 the double occupation number (n^n^) as a function of Coulomb 
interaction energy U at half-filling. It decreases from 1/4 with increasing Coulomb 
interaction, so as to reduce the loss of Coulomb energy U. The MLA-HF state reduces 
more the double occupancy as compared with that of the GW in the weakly correlated 
region, and jumps to the MLA-AA state at U c . In the strongly correlated regime, the 
MLA-AA gives finite value of double occupancy, while the GW gives the Brinkman-Rice 
atom. This verifies the improvement of the Brinkman-Rice atom. The double occupancy 
for the MLA scheme at U c is 0.032, and is consistent with the result of the Quantum 
Monte Carlo (« 0.024), 28 ) though the latter uses the semi-elliptical density of states. 

The momentum distribution for the MLA shows a clear momentum-dependence as a 
function of the HF one electron energy €k a as shown in Fig. 3. It decreases monotonically 
with increasing and shows a jump at the Fermi energy. On the other hand, the 
distribution for the GW is constant below and above the Fermi level. 1_3 ^ The jump 
decreases with increasing U, and disappears beyond U c , indicating the insulating state. 
The curve becomes flatter with further increase of U . 

In summary, we have proposed a new momentum-dependent local ansatz wavefunc- 
tion (MLA) which allows us to describe electron correlations starting from both the 
Hartree-Fock (HF) and the alloy-analogy (AA) limits. The former (i.e., MLA-HF) de- 
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scribes the Fermi-liquid state, while the latter (i.e., MLA-AA) describes the insulator 
state. We have performed the numerical calculations for the half-filled band Hubbard 
model on the hypercubic lattice in infinite dimensions, and demonstrated that the 
ground state energy for the MLA is lower than the GW in the whole range of Coulomb 
interaction. The MLA yields the metal-insulator transition at U c = 3.26. The double 
occupation number is suppressed in the weak and intermediate Coulomb interaction 
regimes as compared with the GW, jumps at U c , and remains finite in the strongly cor- 
related regime as it should be. Finally, we found the momentum distribution functions 
showing a distinct momentum dependence in both the metallic and insulator regimes. 
These results indicate that the MLA approach can overcome the limitations of the 
original MLA, and goes beyond the GW in both the weak and strong U regimes. 

The present work is supported by a Grant-in-Aid for Scientific Research (22540395). 
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